Photonic Simulation of System-Environment Interaction: Non-Markovian Process and 

Dynamical Decoupling 
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The system-environment interaction is simulated by light propagating in coupled photonic waveg- 
uides. The profile of the electromagnetic field provides intuitive physical insight to study the Marko- 
vian and non-Markovian dynamics of open quantum systems. The transition from non-Markovian 
to Markovian process is demonstrated by increasing the size of environment, as the energy evolution 
changes from oscillating to an exponential decay, and the revival period increases. Moreover, the 
dynamical decoupling with a sequence of phase modulations is introduced to such a photonic open 
system to form a band structure in time dimension, where the energy dissipation can be significantly 
accelerated or inhibited. It opens the possibility to tune the dissipation in photonic system, similar 
to the dynamic decoupling of spins. 



Introduction.- Quantum systems always inevitably in- 
teract with environments, which can change the quan- 
tum states and lead to energy dissipation and decoher- 
ence of the systems [l|, Q . For a deep understanding of 
quantum effects and broad apphcation of quantum mat- 
ters, great efforts have been made to study the system- 
environment interaction. Lots of interesting phenom- 
ena of open quantum system have been studied, such 
as the Zeno and anti-Zeno effects [sl, transition be- 
tween quantum Markovian and non-Markovian processes 
[sl, Q . Advanced quantum techniques related to the en- 
vironment have been developed and remarkable progress 
has been achieved in experiments, such as the dissipation 
engineering to prepare and control quantum states with 
the help of environment (7|-[iq|, and dynamical decou- 
pling to preserve the quantum coherence from the envi- 
ronment noise 11-13|. However, it is still difficult to fully 
control the high-dimensional environment, which highly 
limits the understanding, control and application of this 
system-environment interaction. 

Recently, increasing experimental and theoretical ef- 
forts are focused on quantum simulations jlj 
was inspired by Feynman's seminal idea [ij 



|l6|, which 
I . Various 

complex and important physical phenomena can be stud- 
ied with quantum simulators with high efficiency, such as 
the quantum decoherence [l8|, many-body physics, the 
mechanism of superconducting and the general relativ- 
ity [TBI . Such quantum simulations can provide different 
views to study subtle physical processes and reveal new 
phenomena. And more important, people can also learn 
new ideas from the complementary inter discipline, and 
exploit the quantum physics to design new devices for 
quantum technology and practical applications. For ex- 
ample, the photonic simulation of electron spins can be 
used for the topological protected delay [l9| , and the sim- 
ulation of adiabatic passage in waveguide can be applied 
to high-efficiency optical coupler |2QL [2l|. 

In this paper, we simulated the quantum open sys- 
tem by photons in integrated photonic chip. The un- 
derlying physics of the system-environment interaction 



are revealed intuitively by simply observing the electro- 
magnetic field profile. The non-Markovian and Marko- 
vian processes of this photonic open system are studied 
by varying the size of environment. Due to the mem- 
ory effect of finite environment as the energy leaked to 
environment will be reflected back by the boundary of 
environment, the dynamics of the system is non-Markov 
and shows non-exponential decay and revival. While for 
the reservoir with infinite size, the system shows Marko- 
vian exponential decay. Furthermore, we demonstrated 
the control of the system decay by dynamical decoupling. 
With an external modulation of the phase of photon in 
system, the dissipation to environment can be accelerated 
or inhibited. Our study provides an intuitive understand 
of the system-environment interaction, and the results 
can also be used to analyze and reduce the leakage losses 
of photonic structures. 

Model.- As illustrated by the Inset of Fig. 1(a), the 
proposed photonic simulation of the system-environment 
interaction is composed of two separated waveguides. 
Photon can be loaded at system waveguide and travel 
along it (z-axis). The system is open because its energy 
can dissipate to the nearby waveguide which acts as en- 
vironment. The dynamics of photon in single isolated 
waveguide can be described by the Helmholtz equations 
as [22.] 

v^V^(7^) + y(7^)v^(7^) = e^Ij{~^), (1) 

where V{~^) = [e{~^) - 1]^:^ and E = k^, with the di- 
electric relative permittivity e{l^) and the wave- number 
k = 27r/A. The waveguides are uniform along the z-axis. 
Therefore the i-th propagating eigenmode's wave func- 
tion can be expressed as = (^^('^)e*^^^^, where 
(Pi{lt) is the field distribution at the cross section of a 
waveguide, and rii is the effective mode index which can 
be solved from the characteristic equation. For coupled 
system and environment waveguides, we can decompose 
any field distribution perturbatively as 
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FIG. 1: (Color online) (a) The normalized coupling strength 
g^We between system and environment against effect mode 
index for different environment size We. Inset: schematic 
illustration of the coupled photonic waveguides simulation to 
the system-environment interaction, (b) Typical dynamics of 
system energy against the interaction length z for different We , 
with the solid and dashed lines are analytical and numerical 
results, (c)-(e) False-color profiles of electric field intensity 
for light loaded in the system waveguide at z = 0, with d = 
0.15/im and We = 0.23, 5, 10 /xm, respectively. 



where superscript s and e denote system and environ- 
ment waveguides, respectively. Here, i = l,...,7Vs, j = 
1, A^e5 with A^s(e) is the number of modes in the system 
(environment) waveguides. V^|(7^) and ijj^i/^) are eigen- 
modes of waveguides, and c|(z) and c|(z) are correspond- 
ing coefficients. Substituting V{/^) = V^{1^) + V^(^) 
and V^(7^) to Eq. (1), we can obtain the dynamics of the 
modes in each waveguide as 



dz 



j 



)kz eis) 



(3) 



where coupling efficiencies m' 

(e) _ J(^-(-))*y-(-)^;(-)dx 



2n 



and g.) 



2n 



, with all wave functions 



are normalized by / (p^^^"^ dx = 1. 

Considering a thin system waveguide that supports 
only one guiding mode {Ng = 1), while the number of 
modes of the environment waveguide > I which de- 
pends on its width We. Denoting the system and environ- 
ment fields by a = cf e~*^i^^ and bj = glj / gijC^je~^^''3^^ 
{j = 1, A^e) respectively, and replacing the space coor- 



dinate z by time t, we can obtain the Hamiltonian which 
governs the dynamics of this photonic simulation of the 
system-environment interaction as {h = 1) 

H = poa^a -^^Pjb]bj -^^gjia^bj + ab]), (4) 



is the 



where /3o = (nf + m\i)k = n^k and j3j = (n| 
rijk are propagation constants, and gj = y^gfjgfj 
coupling coefficients. Since the direct coupling between 
environment modes is much smaller than other terms, 
mf j {i j) is neglected in following studies. 

Non-Markovian and Markovian Processes.- The 
Hamiltonian of Eq. @ resembles an open system that a 
harmonic oscillator (a) couples to a collection of oscilla- 
tors in environment {bj). Here, both the coupling coeffi- 
cients and the environment can be well controlled, which 
is very promising to simulate the system-environment in- 
teraction. With two-dimensional approximation, Pj and 
gj can be solved analytically. In the following, we stud- 
ied the model with photonic waveguides made of sili- 
con {ud = 3.5), which have been extensively studied in 
practical experiments. The working wavelength (A) is 
1550nm, and the width of system waveguide is fixed to 
Wg = 0.23/im in the single mode regime. 

The normalized coupling strength (g'^We) between sys- 
tem and environment against mode index = Pj/k is 
shown in Fig. 1(a). The continuum of finite size environ- 
ment is quantized to discrete modes, whose density and 
number increases with Wq. The behaviors of g'^We 
converges to a line with infinite A^e when i(;e=oo, which 
corresponds to system-reservoir interaction. With these 
parameters, the dynamics of the system and environment 
can be solved according to Eq. (|4|). In Fig. 1(b), the 
evolutions of energy in system waveguide against z for 
different We are calculated, with light loaded in the sys- 
tem waveguide at z = 0. For a comparison, the fields in 
system are also simulated numerically by finite element 
method. The results of our analytical method agree very 
well with numerical results, with small discrepancies orig- 
inating from the slow varying and weak coupling approx- 
imations used in our analytical expressions. For small 
environment {we = 0.23/im and bfim), the dynamics of 
the system shows periodic oscillation. When the environ- 
ment size increases, the dynamics changes: the system 
energy shows exponential decay at first and revives after 
a distance. 

A direct view of the dynamics of the coupled system- 
environment can be clearly observed in Fig. 1(c)- (e). 
When We is comparable with A, there is few modes can 
be involved to interaction with the system, the coherent 
coupling of which leads to sinusoidal oscillation of sys- 
tem energy [Fig. 1(c)]. For larger more modes could 
interact with the system, giving rise to a complex dynam- 
ics. It is interesting that the corporation of environment 
modes shows a classical ray trajectory of light leaked from 
the system waveguide, as indicated by the dashed line in 
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FIG. 2: (Color online) (a) The memory kernel function for 
different environment sizes, with d = 0.15/im. (b) The dy- 
namics of system for different We and d. (c) The revival period 
R against We, with dots obtained from the system dynamics, 
and solid line from the analytical model, with d — 0.15/im. (d) 
The decay length against d, with dots and line being results 
of the system dynamics and analytical model with We = oo. 
Inset: the decay length against We obtained from the system 
dynamics. 



Fig. 1(d), with an incident angle x = arcsin(/3o//c) with 
respect to z-axis. For We > 10/im, the energy exponen- 
tially decays with z and can almost completely leak into 
environment [Fig. 1(e)]. And after a certain distance, the 
leakage will be reflected into system and cause a revival 
of system. 

Furthermore, some important interaction properties 
to understand the underlying physics can be eas- 
ily obtained from those figures. The formal solu- 
tion of the dynamics of system a{z) is derived in the 
interaction picture Hi{z) = hgj{a^bje'^''^^~^^^^ + 
^j^t^-i{(3o-(3j)z^ by employing the formal solution bj{z) = 
—igj a{r)e~'^^^^~^^'^^dr . The dynamics of system 
reads 



d_ 

dz 



a{z) 



/ a 

Jo 



{t)K{t - z)dr, 



(5) 



with the memory kernel function K{r) = 

small We, the beating of 
multiple environment modes gives rise to the os- 
cillation behavior of K{r). The non-zero memory 
effect leads to non-Markovian dynamics of system, 
showing non-exponential energy decay and revival 
phenomenon [as shown in Fig. 1(b)]. Since the period 
of K{t) is determined by the environment mode density 
p{n) ^ ^^—^ — =, the non-Markovian revival period 

linearly depends on We, which can be observed in 
Fig. 2(b). We can also deduce the revival period 
R = 2wetc[iix + ^0 with intuitive understanding of 
revival as beam reflection at boundary of environment. 



The revival period extracted from the dynamics of 
system shown in Fig. 2(c) are well fitted by this formula. 
The constant Rq is corresponding to an extra distance 
required for tunneling [23| and Goos-Hanchen shift Q 
of light when reflecting at the environment boundary. 

Fig. 2(a) shows that K{t) becomes the Dirac delta 
function 5{t) when the environment size approaches to 
infinity. As a great many environment modes are in- 
volved when We ^ oo, the memory kernel function 
can be approximately written in the form of integration 
K{t) = / g'^{n)e~^^^^~^^^^ p{n)dn^ where the coupling 
strength g{n) is a function of mode index n [Fig. 1(a)]. 
Then 



K{t) ^ J{no)S{r), 
where no = /3o/k and J(n) = g'^{n)we 



(6) 



is the 



spectrum density function of the reservoir. This 5{t) 
memory function gives rise to the Markovian process that 
system energy decays exponentially as a{z) = e~^l^^ . 
The decay length 



L = 



1 



2J(no) 



(7) 



where Lq = ^ 8(ng-i)(n^-ng)3/2 IS solved explicitly. 

In Fig. 2(d), the analytical decay length is compared to 
the results extracted from the dynamics of system with 
logarithmic fitting. Our prediction is consistent with the 
data, and the slight discrepancy is due to non-fiat noise 
spectrum density [Eq. (j6])]. It is surprising that the for- 
mula of L in the continuum limit works for small We, as 
shown in the inset of Fig. 2(d). The oscillation of which 
is due to the variation of mode index of discrete modes 
in finite environment. 

Dynamical Decoupling.- It is well known that external 
modulation to open quantum systems can modify the 
energy decay and decoherence, and the so-called "dy- 
namical decoupling" technique has been widely adopted 



to keep the coherence of electron spins [lll-ll3|. Simi- 
larly, for the photonic simulation of system-environment 
interaction, a dynamical modulation of system waveg- 
uide can also alter the energy decay. Here, a sequence of 
TV modulations with equal interval is applied to the sys- 
tem waveguide, where each modulation corresponds to 
an abrupt change of phase (j). Fig. 3(a) shows the evolu- 
tion of system changes when modulation is applied. The 
dissipation of system can be enhanced or inhibited sig- 
nificantly depending on cj) and N. As shown by Fig. 3(c) 
and (d) are energy in the system waveguide at z = 50/im 
against the modulation phase (/), with different We and 
N. 

There's an intuitive way to understand the modified 
decay: The modulations add extra phase to the propagat- 
ing light, which is equivalent to an increase or decrease of 
the effective index n, as the phase of propagating photon 
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FIG. 3: (Color online) The system evolution (a) and corre- 
sponding spectrum (b) when a sequence of N modulation of 
phase (/) applied to the system, (c) and (d) The system en- 
ergy at 2; = 50iiim against for different We and with 
d — 0.15/im. 



is proportional to n. As shown in Fig. 3(b), the spectrum 
of the system is shifted by the modulations. According to 
the spectrum density of reservoir [Fig. 1(a)], different n 
corresponds to different coupling strength, then the mod- 
ulations give rise to acceleration or deceleration of dissi- 
pation. When n is larger than the cut-off index the 
dissipation of energy in system to the environment is for- 
bidden. From another point of view, the waveguide with 
periodic modulation is similar to the photonic crystal or 
grating structures, which will induce a band structure to 
the light. By this means, noise in environment can be 
prevented from propagating in system. These provide 
physics insight to the dynamical decoupling in time di- 
mension where a sequence of modulation in time-axis is 
applied: the effective frequency of the system shifts to a 
higher value than the cut-off frequency of bath. Or, we 
can image a crystal or band structure in time dimension, 
where only the signal with certain frequencies can enter 
and be kept in the system while most of the broadband 
noises in environment is blocked. 

One possible way to realize such phase modulation se- 
quence is shown in Fig. 4(a), with whispering gallery 
mode (WGM) microdisks coupling to the system waveg- 
uide [251 . The microdisk with radius r = Ijim has intrin- 
sic quality factor higher than 10^. When put close to the 
system waveguide with a gap (0.15/im), its loaded quality 
factor is only about 2000. That means the WGM is work- 
ing in the strongly over coupling regime with ^ a^^, 
where /^^(e) is the intrinsic (external) loss. When the light 
in system passes the over-coupled WGM resonator, the 
change of transmitted light is T(A) = — ^ 

— 1^-^/^'' , where A is the frequency detuning to the reso- 
nance. As |T(A)| ^ 1, the system acquires a modulation 
of phase ^(A) = arg[T(A)]. Numerical simulation of 
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FIG. 4: (Color online) (a) The fields of system and environ- 
ment when microdisks are added to the waveguide to induce 
a phase modulation, (b) The dynamics against the detuning 
of input light to the whispering gallery modes. 



the dynamical decoupling with = 10 for = 10/im 
is performed. Comparing the mode profile of Fig. 4(a) 
with Fig. 1(e), the evolutions are significantly changed 
by modulations. From Fig. 4(b), a strong modification 
of the decay is shown around the resonance that phase 
modulation ^ tt, which is consistent with the predic- 
tion in Fig. 3(a). 

Conclusion.- The Markovian, non-Mar kovian pro- 
cesses and dynamic decoupling of open quantum sys- 
tems are studied in a photonic simulation of system- 
environment interaction. Intuitive physical insight and 
deep understanding of these phenomena can be gained 
from the direct view of electromagnetic field profiles. Our 
study also opens the possibility to tune the dissipation 
in photonic system, similar to the dynamic decoupling of 
spins. 
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